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Current  System  (LCS) ,  several  experiments  are  conducted 
with  a  sigma  coordinate  primitive  equation  model  on  a  beta- 
plane.  Results  show  that  the  LCS  is  an  anomalous  eastern 
boundary  current  (EEC)  that  generates  a  coastal  poleward 
current,  an  equatorward  undercurrent,  and  highly  energetic 
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intensifying  and  trapping  currents  near  the  coast, 
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I. 


INTRODUCTION 


Subtropical  eastern  boundary  currents  (EEC)  are 
usually  equatorward  at  the  surface,  with  poleward 
undercurrents,  which  are  forced  by  prevailing  equatorward 
winds.  The  major  EBCs  such  as  Peru,  California,  Benguela 
and  Canary  are  distinguished  by  surface  dynamic  height 
fields  which  decrease  toward  the  coast  and  near-surface 
isopycnals  that  slope  upward  (Wooster  and  Reid,  1963; 
McCreary  et  al .  ,  1986).  These  systems  are  characterized  by 
climatologically  weak  (<10  cm/s),  broad  (-1000  km  wide) 
surface  flow  toward  the  equator,  cold  upwelled  water  at  the 
surface,  shallow  (<30  m  depth)  thermoclines ,  and  high 
biological  production  (Parrish  et  al . ,  1983)  due  to  regions 
of  significant  upwelling.  They  are  one  component  of  the 
subtropical  anticyclonic  gyres,  which  are  driven  primarily 
by  anticyclonic  wind  fields. 

Along  the  coast  of  Western  Australia,  the  prevailing 
winds  are  predominantly  equatorward  (Thompson,  1984; 
Godfrey  and  Ridgeway,  1985);  however,  unique  to  the  region, 
is  a  poleward  surface  current  known  as  the  Leeuwin  Current. 
Observational  studies  along  the  western  coast  of  Australia 
have  shown  that  this  current  is  characterized  by  a  strong 
(>150  cm/s  at  times),  narrow  (<100  km  wide),  poleward 
surface  current  that  flows  opposite  the  prevailing  wind 
direction  (Cresswell  and  Golding,  1980;  Godfrey  et  al .  , 
1986),  anomalous  warm  water  at  the  surface,  a  deep  (>50  m 
depth)  thermocline  (Thompson,  1984)  and  lower  biological 
production  due  to  vast  regions  of  downwelling  (Batteen  et 
al . ,  1992) . 
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An  anomalous  EBC,  the  Leeuwin  Current  is  driven  by  the 
uncharacteristically  large  thermohaline  gradient  along  the 
western  Australian  coast.  The  current  originates  near 
Shark  Bay  and  flows  poleward  along  the  continental  shelf 
(~200  m  depth)  off  the  coast  of  western  Australia  to  Cape 
Leeuwin  (see  Figure  1  for  geographical  locations)  and  then 
eastward  into  the  Great  Australian  Bight  (Cresswell  and 
Golding,  1980)  .  Associated  with  the  system  are  near¬ 
surface  isopycnals  that  slope  downward,  surface  dynamic 
height  fields  that  increase  toward  the  coast,  and  an 
equatorward  undercurrent  (Thompson,  1984;  McCreary  et  al .  , 
1986)  .  The  equatorward  wind  stress  is  overwhelmed  by  the 
meridional  pressure  gradient  created  from  excessive  heating 
in  the  equatorial  region  and  large  amounts  of  cooling  in 
the  poleward  region  (Godfrey  and  Ridgway,  1985),  generating 
the  Leeuwin  Current.  The  source  for  the  Leeuwin  Current  is 
predominantly  an  alongshore  steric  height  gradient  due  to 
tropical  Pacific  water  from  the  Indonesian  throughflow 
(Godfrey  and  Ridgway,  1985;  Hirst  and  Godfrey,  1993), 
augmented  by  geostrophic  inflow  from  the  west  (McCreary  et 
al . ,  1986;  Thompson,  1987)  .  The  poleward  flow  of  the 
Leeuwin  Current  intensifies  due  to  the  strong  inflow  of  the 
subtropical  waters  toward  the  coast  (Batteen  et  al . ,  1992). 

Below  the  Leeuwin  Current,  there  is  an  anomalous 
equatorward  current  off  western  Australia  (Church  et  al . , 
1989)  .  This  undercurrent  is  centered  near  450  m  depth  and 
can  attain  speeds  comparable  to  the  surface  flow,  exceeding 
on  average  10  cm/s  (Smith  et  al .  ,  1991)  .  Smith  et  al . 
(1991)  observed  core  speeds  of  the  narrow  undercurrent  up 
to  30  cm/ s  between  ~250  and  350  m  depth  and  found  that  the 
current  was  usually  confined  to  the  continental  shelf  slope 
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between  ~250  and  450  m  depth.  Although  there  is  evidence 
of  a  westward  flow  of  ~20  cm/ s  centered  between  -400  and 
700  m  depth  off  the  southern  coast  at  116°  E  (Cresswell  and 
Peterson,  1993),  no  studies  to  date  have  clearly  resolved 
whether  there  is  a  westward  undercurrent  off  Southern 
Australia  . 

Previous  numerical  modeling  studies  by  Batteen  et  al . 
(1992)  investigated  the  effects  of  annual  climatological 
wind  forcing  and  initialized  thermohaline  gradients  on  the 
Leeuwin  Current  System  (LCS) ,  but  the  study  was  limited  to 
the  coast  off  western  Australia  and  did  not  include  the 
influence  of  topography.  Batteen  and  Huang  (1998)  studied 
the  effects  of  salinity  on  the  density  of  the  LCS.  Batteen 
and  Butler  (1998)  examined  the  effects  of  continuously 
forced  annual  Indian  Ocean  thermohaline  gradients  on  the 
LCS  and  extended  the  domain  to  include  both  the  western  and 
southwestern  coasts  of  Australia.  Batteen  et  al .  (2000) 
investigated  the  effects  of  seasonal  forcing  on  the  LCS 
with  a  full  primitive  equation  ocean  model  but  again  only 
considered  flat  bottom  cases. 

This  study  seeks  to  extend  prior  efforts  in  this  area 
by  including  the  addition  of  realistic  topography,  allowing 
a  better  understanding  of  the  LCS  by  including  the 
influences  of  topography.  The  Princeton  Ocean  Model  (POM), 
a  bottom  following,  sigma  coordinate  model,  was  chosen  for 
this  study  because  it  has  been  widely  used  to  simulate 
coastal  processes  associated  with  continental  shelf  flows 
and  bottom  boundary  layer  dynamics.  The  results  of  several 
numerical  experiments  (Table  1)  are  explored.  Each 
experiment  is  run  on  a  beta  plane.  In  Experiment  1  the 
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horizontally  averaged  annual  climatology  is  used  with 
annual  wind  forcing  on  a  flat  bottom.  Experiment  2 
investigates  the  effect  of  thermohaline  gradient  forcing  by 
using  full  annual  climatology  and  no  wind  forcing  on  a  flat 
bottom.  The  final  flat  bottom  experiment.  Experiment  3 
includes  both  wind  and  thermohaline  gradient  forcing 
mechanisms  to  gain  an  understanding  of  the  combined  effect 
of  the  forcing  mechanisms  on  the  dynamics  of  the  UTS.  In 
Experiment  4,  velocity  errors  produced  by  the  pressure 

gradient  force  error  (an  error  inherent  in  all  three- 
dimensional,  sigma  coordinate  models)  in  the  presence  of 
bottom  topography  are  determined  using  horizontally 

averaged  annual  climatology  with  bottom  topography  and  no 
wind  forcing.  Experiments  1  through  3  are  repeated  as 
experiments  5,  6  and  7,  respectively,  but  with  the  addition 

of  bottom  topography. 

Note  that  the  results  from  experiments  with  wind 
forcing  only,  experiments  1  and  5,  are  more  typical  of 
classical  EBCs  than  of  the  LCS .  They  are  presented  to 

highlight  to  effect  of  wind  forcing  in  the  absence  of  a 

thermohaline  gradient.  In  contrast  the  final  experiment, 
which  includes  bottom  topography  with  annual  thermohaline 
gradient  and  wind  forcing,  provides  the  most  accurate 
representation  of  the  LCS. 


This 

study  is  organized  as 

follows . 

In 

section  II  the 

numerical 

model  is  described 

along 

with 

the 

specific 

experimental  conditions.  The 

results 

of 

the 

numerical 

experiments  are  presented  in  section  III,  along  with  a 
discussion  comparing  the  results  to  observations.  A 
summary  of  the  study  is  provided  in  section  IV. 
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II .  MODEL  DESCRIPTION 


A .  DATA  SETS 

The  model  domain  (Figure  1)  encompasses  the  western 
and  southwestern  coasts  of  Australia,  from  24.5°S  to  39°S, 
and  from  109°W  to  121°W.  The  topographic  data  was  obtained 
from  the  Institute  of  Geophysics  and  Planetary  Physics, 
University  of  California  San  Diego  (Sandwell  and  Smith, 
1996)  .  The  bottom  topography  data  set,  compiled  from  over 
30-years  of  bottom  echo  sounding  by  ships,  has  a  resolution 
of  2  minutes  (1/30  of  a  degree)  .  Altimetry  data  was  used 
to  interpolate  soundings  in  data  sparse  regions. 

Annual  temperature  and  salinity  values  were  obtained 
from  Levitus  and  Boyer  (1994)  and  Levitus  et  al .  (1994)  . 
These  data  sets  incorporate  a  1  by  1  degree  horizontal 
resolution  at  the  thirty-three  vertical  levels  shown  in 
Table  2 . 

B .  PRE-PROCESSING 

The  original  topography  (Figures  2a  and  2b)  was 
interpolated  to  the  resolution  used  in  the  POM  Model,  i.e., 
10  km  by  10  km  offshore,  3  km  by  10  km  alongshore  and  3  km 
by  3  km  around  the  southwest  corner  of  Australia  with  a 
total  of  252  by  226  grid  points  (Figure  3)  .  The  highest 
resolution  was  used  where  the  values  of  the  slope  parameter 

15// 1  _ 

(defined  by  Mellor  et  al . ,  1998,  as  - where  H  is  the 

2*H 

average  depth  and  dH  is  the  difference  in  depth  between  two 
adjacent  cells)  ,  were  the  largest  in  both  the  latitude  and 
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longitude  directions.  The  higher  resolution  bands  near  the 
southwest  corner  of  Australia  were  used  to  minimize  the 

slope  parameter  (Figure  3)  .  Since  over  much  of  the  area 
the  slope  parameter  was  greater  than  0.2,  which  is  the 
suggested  maximum  value  to  be  used  in  a  sigma  coordinate 

models  (Mellor  et  al.,  1998),  the  topography  was  first 

smoothed  with  a  Gaussian  filter.  The  new  depth  of  each 
point  calculated  with  this  filter  is  a  Gaussian-shaped, 

weighted  average  of  25  by  25  points  with  a  standard 
deviation  of  8.  Subsequent  depths  greater  than  2500  m  were 
reassigned  to  depths  of  2500m,  with  land  assigned  a  depth 
of  20  m  (to  avoid  divisions  by  zero  in  the  model)  .  The 

resulting  topography  is  shown  in  Figures  4a  and  4b. 

The  annual  temperature  and  salinity  values  were 

interpolated  for  the  horizontal  spatial  resolution  of  the 
model  and  for  the  21  vertical  sigma  levels  (Table  3  and 

Figure  5)  with  a  three-dimensional  (3D)  linear 

interpolation  scheme.  This  had  to  be  done  separately  for 
smoothed  topography  and  for  flat  bottom  topography  due  to 
the  change  in  vertical  levels  between  the  two  data  sets. 
Plots  of  annual  temperature  and  salinity  fields  for  the 
surface  (sigma  level  one)  are  shown  in  Figures  6  and  7, 
respectively . 

To  obtain  an  annual  non-weighted  average  wind  vector 
field  (Figure  8),  ten  years  of  daily  seasonal  winds  were 
averaged.  To  be  compatible  with  the  other  data  fields,  the 
wind  vectors  were  interpolated  to  the  same  horizontal 
spatial  resolution  of  the  model  with  a  2D  linear 
interpolation  scheme. 
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C.  BRIEF  MODEL  DESCRIPTION 

The  POM,  a  well  documented  model  (e.g.,  Blumberg  and 
Mellon,  1987;  Mellon,  1996),  was  used  in  the  model  studies. 
POM  is  a  pnimitive  equation,  fnee  sunface  model  with  a 
second-moment  tunbulence  closune  scheme  (Mellon  and  Yamada, 
1982)  that,  thnough  the  use  of  bottom-following  sigma 
levels,  can  nealist ically  simulate  pnocesses  associated 
with  continental  shelf  flows  and  bottom  boundany  layen 
dynamics  in  local  domains  (e.g.,  bays,  estuanies  and 
coastal  negions) .  Recently,  the  model  has  been  used 
successfully  to  simulate  decadal  pnocesses  in  entine  ocean 
basins  (see  Ezen  and  Mellon,  1994,  1997)  . 

As  descnibed  eanlien,  the  nesolution  of  the  honizontal 
onthogonal  gnid  vanies  between  3  km  by  3  km  and  10  km  by  10 
km  (Figune  3)  .  The  vaniable  gnid  allows  the  use  of  mone 
(less)  points  in  negions  of  lange  (small)  gnadients. 

The  21  ventical  sigma  levels  used  ane  shown  in  Figune 
5  and  Table  3.  The  sigma  values  nange  fnom  zeno  at  the 
sunface  to  minus  one  at  the  bottom  with  the  ventical  gnid 
spacing  pnopontional  to  the  ocean  depth.  The  ventical 
nesolution  has  been  chosen  to  be  highen  nean  the  sunface 
and  the  bottom  in  onden  to  nesolve  both  the  sunface 
boundany  layen  and  the  bottom  boundany  layen,  which  ane 
impontant  in  coastal  negions.  To  eliminate  the  time 
constnaints  fon  the  ventical  gnid  nelated  to  highen 
nesolution  nean  sunface,  bottom  and  shallow  waten,  an 
implicit  ventical  time  diffenencing  scheme  is  used. 

The  pnognostic  vaniables  of  the  model  ane  potential 
tempenatune,  salinity,  density,  the  thnee  velocity 
components,  sunface  elevation,  tunbulent  kinetic  enengy  and 
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length  scale.  The  model  has  a  split  time  step  interval  for 
external  and  internal  modes.  The  external  mode  solves  the 
equations  for  the  vertically  integrated  momentum  equations. 
It  also  provides  the  sea  surface  and  barotropic  velocity 
components,  and  has  a  time  step  of  6  seconds.  The  internal 
mode  solves  the  complete  3D  equations  and  has  a  time  step 
of  300  seconds. 

A  Smagorinsky  formulation  (Smagorinsky  et  al . ,  1965) 

is  used  for  the  horizontal  diffusion  in  which  the 
horizontal  viscosity  coefficients  depend  on  the  grid  size, 
the  velocity  gradient  and  a  coefficient.  In  this  study  a 
value  of  0.2  was  assigned  to  this  coefficient,  consistent 
with  other  POM  studies  (e.g.,  Ezer  and  Mellor,  1997) . 

D.  INITIALIZATION,  FORCING  AND  BOUNDARY  CONDITIONS 

The  model  was  initialized  with  annual  temperature  and 
salinity  values  obtained  from  Levitus  and  Boyer  (1994)  and 
Levitus  et  al .  (1994)  .  Since  the  model  runs  reached  a 

quasi-equilibrium  state  in  a  relatively  short  time  (~60 
days),  zero  salinity  and  temperature  fluxes  were  prescribed 
at  the  ocean  surface.  The  annual  climatological  surface 
temperature  (Figure  6)  shows  a  strong  north-south  gradient 
with  a  maximum  temperature  near  Shark  Bay  of  ~24°C  and  a 
minimum  of  14°C  at  the  southern  edge  of  the  model  domain. 
The  climatological  surface  salinity  (Figure  7)  shows  a 
maximum  of  ~38.9  at  the  center  of  the  western  edge  of  the 
model  domain  and  a  minimum  of  ~35.3  along  the  northern 
Australian  coast.  Models  using  horizontally  averaged 

climatology  use  the  same  initial  temperature  and  salinity 
over  the  entire  model  domain  (e.g,  19.8°C  and  35.39  are  the 


temperature  and  salinity  values,  respectively,  for  the 
surface) . 

The  model  was  forced  from  rest  with  the  annual 
European  Center  for  Medium-range  Weather  Forecasting 
(ECMWF)  wind  fields,  which  were  interpolated  for  the  model 
grid.  As  seen  in  Figure  8,  westerly  winds  dominate  the 
southern  portion  of  the  model  domain.  Along  the  west  coast 
of  Australia  the  wind  are  southerly  to  south-southeasterly 
and  generally  increase  in  strength  away  from  the  coast. 

Correct  specification  of  the  open  boundary  conditions 
(BC)  is  important  to  achieve  realistic  results,  with  no 
reflections,  clamping,  spurious  currents  or  numerical 
alteration  of  the  total  volume  of  water  in  the  model.  The 
problem  is  that  there  is  not  a  general  criterion  that  can 
give  the  answer  to  what  boundary  conditions  are  the  best 
for  a  specific  model  or  study.  For  models  with  a  free 
surface,  such  as  used  here,  one  of  the  important  criteria 
is  that  the  BCs  should  be  transparent  to  the  waves.  In 
this  model,  a  gradient  boundary  condition  (Chapman,  1985), 
which  allows  geostrophic  flow  normal  to  the  boundary, 
worked  best  for  the  elevation.  For  baroclinic  velocity 
components  normal  to  the  boundary,  an  explicit  wave 
radiation  scheme  based  on  the  Sommerfield  radiation 
conditions  was  used.  For  inflow  situations,  the  model  was 
forced  with  annual  temperature  and  salinity  values  obtained 
from  Levitus  and  Boyer  (1994)  and  Levitus  et  al .  (1994), 
while  in  outflow  situations  an  advection  scheme  was  used. 
An  improvement  on  previous  models,  is  the  introduction  of 
an  advanced  volume  constraint  subroutine  based  on  research 
performed  by  Marchesiello  et  al . ,  2001,  which  drastically 
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reduced  the  loss  of  volume  at  the  boundaries  to 
insignificant  levels. 

For  the  barotropic  velocity  components,  a  Flatter 
radiation  plus  Roed  local  solution  (FRO)  was  used.  Palma 
and  Matano  (2000)  showed  good  results  with  the  FRO  solution 
during  BC  tests  to  determine  the  BCs  response  to  an 
alongshelf  wind  stress.  Palma  and  Matano  (1998)  also 
showed  that  the  FRO  BC  demonstrated  good  reflection 
properties  and  results  in  a  test  that  determined  the  BC 
response  to  the  combined  action  of  wind  forcing  and  wave 
radiation.  Their  tests  were  executed  with  the  barotropic 
version  of  POM  and  compared  with  the  benchmark  results  (no 
boundary  conditions) . 
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III.  RESULTS  FROM  MODEL  SIMULATIONS 


A.  WIND  FORCING  ON  A  FLAT  BOTTOM 

In  Experiment  1  (see  Table  1),  the  model  was 
initialized  with  the  horizontally  averaged  annual 
climatological  temperature  and  salinity.  A  realistic 
coastline  and  flat  bottom  (constant  depth  of  2500  m)  were 
used,  and  the  model  was  forced  with  annual  climatological 
winds.  The  goal  of  this  experiment  was  to  highlight  the 
role  of  wind  forcing  alone  in  the  LCS . 

As  expected,  the  southerly  winds  along  the  coast  have 
generated  a  weak  equatorward  current  of  ~20  cm/s  (Figure 
9a)  .  By  day  20,  the  equatorward  winds  induce  net  offshore 
surface  Ekman  transport,  which  generates  regions  of 
upwelling  evidenced  by  the  cool  water  along  the  coast  near 
Shark  Bay.  The  most  intense  upwelling  is  associated  with 
the  stronger  winds  in  the  north,  but  the  upwelling  effect 
is  apparent  all  along  the  western  coast  of  Australia.  A 
typical  cross-section  of  meridional  velocity  (Figure  9b)  , 
shows  the  core  of  the  poleward  undercurrent  at  ~600  m  depth 
with  a  speed  of  ~2  cm/s  underlying  the  equatorward  surface 
current  with  a  core  speed  of  ~16  cm/ s . 

By  day  40  (Figure  9c) ,  the  coastal  upwelling  has 
intensified  bringing  18°C  water  to  the  surface  near  Shark 
Bay.  In  addition,  the  offshore  extent  of  the  coastal 
upwelling  is  much  wider  than  at  day  20.  As  expected,  capes 
and  promontories  enhance  the  upwelling  regions  to  the 
north.  A  cold  core,  cyclonic  eddy  has  formed  near  Shark 
Bay.  The  velocity  of  the  equatorward  surface  current 
remains  at  ~20  cm/s,  implying  that  the  wind  has  already 
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imparted  its  maximum  amount  of  energy  into  the  current. 
The  width  of  the  surface  current  has  increased  to  ~50  km. 
Simultaneously,  as  the  undercurrent  intensifies  to  ~5  cm/s, 
it  shallows  and  displaces  the  equatorward  current  offshore 
(e.g.,  see  Figure  9d) . 

By  day  60  (Figure  9e)  ,  the  effects  of  upwelling  are 
evident  along  the  entire  western  coast  and  extend  south 
around  Cape  Leeuwin.  Meanders,  eddies  and  filaments  are 
much  more  apparent  in  the  northern  end  of  the  domain  where 
the  upwelling  is  the  most  intense.  Inside  Shark  Bay  the 
water  temperature  has  decreased  about  2°C  from  the  initial 
value  at  the  start  of  the  model  run.  The  north-south 
velocity  profile  (Figure  9f ) ,  shows  an  eddy  forming  at 
34.3°S,  evidenced  by  the  reversal  in  surface  current 
direction . 

The  continued  widening  of  the  equatorward  coastal  jet 
(e.g.,  see  Figures  9b,  9d  and  9f)  is  consistent  with  the 
results  of  McCreary  et  al .  (1987)  and  Batteen  et  al. 
(1989) .  They  showed  that,  due  to  the  beta  (h )  effect,  the 
surface  coastal  jet  does  not  necessarily  have  to  be 
confined  to  within  a  Rossby  radius  of  deformation  of  the 
coast.  The  h  effect  also  allows  the  existence  of  freely 
propagating  planetary  waves,  i.e.,  Rossby  waves  (Gill, 
1982)  .  The  offshore  propagation  of  these  waves  contributes 
to  the  generation  of  an  alongshore  pressure  gradient  field, 
which  can  aid  the  development  of  subsurface  currents  along 
the  eastern  boundary.  As  a  result,  the  h  effect  changes 
both  the  vertical  and  horizontal  structure  of  the  surface 
and  subsurface  currents.  As  the  undercurrent  intensifies 
and  shoals,  strong  vertical  and  horizontal  shear  occurs  in 
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the  upper  layers  creating  baroclinically  and  barotropically 
unstable  coastal  currents.  This  unstable  condition  leads 
to  the  development  of  the  meanders,  eddies  and  filaments 
evident  in  the  northern  end  of  the  model  domain. 

The  results  from  this  experiment  resemble  the  dynamic 
characteristics  associated  with  classical  EBCs,  where  wind 
forcing  is  the  predominant  forcing  mechanism.  Upwelling 
was  produced  due  to  favorable  winds;  however,  these  results 
were  not  representative  of  the  LCS .  Additional  experiments 
will  study  other  factors  that  distinguish  the  LCS  from 
other  EBC  regions. 

B.  THERMOHALINE  GRADIENT  FORCING  ON  A  FLAT  BOTTOM 

Experiment  2  (see  Table  1)  focuses  on  the  thermohaline 
gradient.  The  experiment  is  initialized  with  annual 
climatology  that  has  not  been  horizontally  averaged  but 
without  wind  forcing.  It  uses  the  same  realistic  coastline 
and  flat  bottom  as  in  Experiment  1.  This  experiment 
introduces  the  conditions  unique  to  the  LCS. 

In  stark  contrast  to  Experiment  1,  by  day  20  (Eigure 
10a)  an  energetic  and  dynamic  poleward  current  has 
developed  with  a  maximum  velocity  of  ~120  cm/s .  The 
poleward  flow  is  evident  along  the  entire  western  coast  of 
Australia  and  extends  around  Cape  Leeuwin  into  the  Great 
Australian  Bight.  As  the  current  rounds  Cape  Leeuwin  it  is 
steered  eastward  due  to  the  Coriolis  effect.  As  the 
current  travels  south,  it  advects  warm,  subtropical  water 
with  temperatures  off  the  coast  of  Cape  Leeuwin  reaching 
21°C.  The  strength  of  the  surface  current  is  created  by  the 
large  north-south  thermohaline  gradient.  The  thermohaline 


13 


gradient  generates  a  strong  onshore  geostrophic  inflow.  As 
the  inflow  nears  the  coast,  it  is  deflected  southward 
creating  a  narrow  core  of  warm  water  along  the  entire  coast 
of  Western  Australia.  Continual  onshore  flow  augments  the 
poleward  flow  throughout  the  model's  duration.  In  addition 
to  a  strong  surface  current,  an  intense  equatorward 
undercurrent  is  established  (e.g.,  see  Figure  10b)  .  The 
undercurrent  is  centered  near  600  m  depth  with  a  maximum 
velocity  of  ~20  cm/s,  and  underlies  the  equatorward  surface 
current  with  a  core  speed  of  ~  80  cm/s . 

By  day  40  (Figure  10c) ,  the  poleward  surface  current 
has  increased  to  ~140  cm/ s .  As  was  seen  in  Experiment  1, 
the  width  of  the  current  increases  due  to  the  P  effect. 
Warm,  subtropical  water  has  been  advected  farther  around 
Cape  Leeuwin,  with  surface  temperatures  off  the  Cape 
reaching  25°C.  The  barotropic  instability  of  the  current  is 
enhanced  by  the  horizontal  temperature  gradient  and  the 
increased  strength  of  the  current  producing  large 
horizontal  shear.  In  addition,  the  undercurrent  (Figure 
lOd)  is  intensifying,  ~25  cm/s,  and  shoaling,  increasing 
the  baroclinic  instability.  Large  meanders  form  south  of 
Perth  and  around  the  coast  to  the  east  due  to  coastline 
irregularities  and  the  instability  of  the  current.  Near 
Cape  Leeuwin  cold,  sub-polar  water  is  being  entrained  west 
of  the  main  flow  generating  dipole  eddy  pairs,  which 
subsequently  propagate  eastward  along  the  southern  coast  of 
Australia . 

By  day  60  (Figure  lOe)  ,  the  dynamics  of  the  current 
resemble  those  at  day  40.  Off  of  Perth  the  current 
continues  to  broaden  due  to  P  effects.  The  entrainment  of 
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cold,  subpolar  water  is  evident  in  the  large  eddy  dipole 
pairs.  The  dipole  pairs  are  relatively  flat  and  broad  due 
to  the  (3  effect  stretching  them  westward  as  they  propagate 
eastward  along  the  southern  coast.  The  current  has 
advected  warm  25°C  water  well  south  of  Cape  Leeuwin.  A 
north-south  cross  section  (Figure  lOf)  dissects  the  eddy 
forming  off  Cape  Leeuwin.  The  eddy  is  centered  near  34.3°S, 
113. 5°W,  with  the  broad,  poleward  surface  current  present 
along  the  coast. 

The  results  of  this  experiment  are  more  characteristic 
of  the  LCS,  generating  a  strong  poleward  surface  current 
with  an  equatorward  undercurrent.  It  demonstrates  the 
importance  of  the  strong  climatological  thermohaline 
gradient.  However,  the  model  currents  were  too  broad  along 
the  coast  and  the  eddies  were  not  very  well  defined. 

C.  WIND  AND  THERMOHALINE  GRADIENT  FORCING  ON  A  FLAT 

BOTTOM 

The  final  flat  bottom  experiment.  Experiment  3  (see 
Table  1),  most  closely  simulates  the  LCS  by  using  both  wind 
and  thermohaline  forcing.  The  model  was  initialized  with 
full  annual  temperature  and  salinity  climatology.  In 
addition,  annual  wind  forcing  was  added  to  determine  the 
net  result  of  both  forcing  mechanisms. 

Similar  to  the  no  wind  case,  by  day  20  (Figure  11a),  a 
warm,  poleward  surface  current  has  formed  along  the  coast; 
however,  because  the  wind  and  thermohaline  forcing  oppose 
each  other,  the  current  is  not  as  strong  as  in  Experiment  2 
and  has  a  maximum  velocity  of  ~90  cm/ s .  Overall  the 
dynamics  and  temperature  advection  of  the  current  are  very 
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similar  to  the  results  found  at  day  20  of  Experiment  2 
(Figure  10a).  Additionally,  a  realistic  equatorward 
undercurrent  has  formed  off  Cape  Leeuwin,  Lat  34.3°S, 
centered  near  500  m  depth  with  a  speed  of  20  cm/s  (Figure 
11b)  . 

By  day  40  (Figure  11c)  ,  the  current  speed  has 
increased  to  120  cm/s,  consistent  with  the  speed  in 
Experiment  2.  Again,  compared  to  the  no  wind  experiment 
(Experiment  2),  the  meanders  and  eddy  formation  regions  are 
similar,  with  the  exception  of  a  decrease  in  the  speed  of 
the  surface  current.  A  typical  cross-section  (Figure  lid) 
also  shows  that  the  current  does  not  have  as  wide  a  coastal 
extent  as  before.  This  can  be  attributed  to  the  addition 
of  wind  forcing  that  opposes  the  thermohaline  forced 
current  and  reduces  its  offshore  extent. 

Day  60  (Figure  lie)  shows  that  although  the  velocity 
has  not  increased,  the  current  broadens  north  of  Perth, 
while  also  continuing  to  meander  and  create  large  eddies 
off  Cape  Leeuwin.  The  addition  of  wind  forcing  helps  to 
drive  equatorward  flow  in  dipole  eddy  pairs  making  them 
more  defined  than  in  Experiment  2  .  A  meridional  cross- 
section  (Figure  Ilf)  shows  generation  of  an  eddy  off  Cape 
Leeuwin.  The  core  of  the  current  has  moved  off  the  coast 
and  defines  the  outer  edge  of  the  eddy,  with  the  poleward 
(blue  shading)  flow  inshore  of  the  main  current.  A 
meridional  cross-section  farther  equatorward  (Figure  llg) 
shows  that  as  the  main  surface  current  leaves  the  coast  off 
Perth  (due  to  meanders  or  coastline  irregularities)  ,  wind 
forcing  dominates  and  creates  an  equatorward  current.  As 
expected,  the  equatorward  flow  generates  regions  of 
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upwelling  which  are  distinguished  by  the  slightly  cooler 
surface  temperatures  off  Perth  and  to  the  north  of 
promontories,  such  as  at  ~28°S. 

D.  PRESSURE  GRADIENT  FORCE  ERROR  DETERMINATION 

In  Experiment  4,  the  model  was  initialized  at  rest 
with  horizontally  averaged  annual  climatological 
temperatures  and  salinities.  A  realistic  coastline  and 
bottom  topography  were  used  but  no  wind  or  thermohaline 
forcing  was  permitted. 

With  horizontally  averaged  climatology  and  no  forcing 
mechanisms,  we  would  expect  that  nothing  will  happen,  i.e., 
the  initial  state  of  rest  should  be  maintained.  Due  to 
pressure  gradient  force  errors  (explained  below) ,  however, 
this  will  not  be  the  case  and  there  will  be  velocities  that 
result  from  these  errors. 

Velocity  errors  induced  by  the  pressure  gradient  force 
are  unavoidable  in  3D  sigma-coordinate  models.  Two  types 
of  sigma-coordinate  errors  exist,  the  sigma  error  of  the 
first  kind  (SEEK)  and  of  the  second  kind  (SESK) ,  as  defined 
by  Mellor  et  al .  (1998)  .  The  first  one  goes  to  zero 
prognost ically  by  advecting  the  density  field  to  a  new 
state  of  equilibrium.  The  second  one,  a  vorticity  error, 
is  the  most  important  because  it  does  not  vanish  with  time, 
and  is  present  in  both  2D  and  3D  cases. 

There  are  several  techniques  to  reduce  the  pressure 
gradient  errors : 

1 .  Smoothing  the  topography  can  reduce  both  SEEK  and 

SESK.  In  particular,  the  slope  parameter  should  not 
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be  greater  than  0.2  (Mellor  et  al .  , 


1998) .  Greater 


values  can  artificially  induce  currents  over  1  m/s. 

2.  Using  the  highest  possible  resolution  can  reduce 
the  errors,  since  the  pressure  gradient  error 
decreases  with  the  square  of  the  horizontal  and 
vertical  grid  size. 

3.  Subtracting  the  horizontally  averaged  density 
before  the  computation  of  the  baroclinic  integral 
reduces  the  SESK. 

4.  Using  a  curvilinear  grid  that  follows  the 
bathymetry  reduces  the  SESK. 

This  study  used  the  first  three  techniques.  The  use 
of  a  curvilinear  grid  was  not  employed  since  the  first 
three  techniques  successfully  reduced  the  error  to  an 
acceptable  level  (to  ~1  cm/s  or  less) . 

To  show  where  the  remaining  velocity  errors  (induced 
by  the  pressure  gradient  force  error)  are  present  in  the 
model  domain,  the  velocity  field  for  day  60  at  sigma  level 
1  (surface)  is  shown  in  Eigure  12.  Maximum  velocities  of 
up  to  ~1  cm/ s  are  found  along  the  coast  where  the  slope 
parameter  is  the  greatest.  These  results  are  similar  to 
those  obtained  by  Martinho  (2001)  and  indicate  that  with 
the  use  of  the  three  error  reduction  techniques,  the 
pressure  gradient  force  error  has  been  considerably  reduced 
from  values  over  ~1  m/s  to  values  less  than  1  cm/s. 
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E.  WIND  FORCING  OVER  TOPOGRAPHY 

Experiment  5  (see  Table  1)  mimics  Experiment  1,  except 
instead  of  a  flat  bottom,  topography  was  added  with  the 
same  model  forcing  mechanisms.  The  model  was  initialized 
with  horizontally  averaged  annual  climatological 
temperature  and  salinity.  The  same  coastline  and  annual 
climatological  wind  field  were  used.  The  results  of 
Experiment  5  were  run  to  isolate  the  effects  of  topographic 
beta  and  demonstrate  its  role  in  a  wind-forced  model. 

By  day  20  (Eigure  13a),  an  equatorward  current  has 
developed  with  a  maximum  speed  of  ~40  cm/ s .  Weak, 
localized  coastal  upwelling  is  evident  near  promontories 
north  of  Perth.  The  upwelling  is  not  as  extensive  or 
widespread  as  in  the  flat-bottom  case  (Experiment  1)  .  A 
typical  cross-section  (Eigure  13b)  shows  that  the  core  of 
the  surface  current  is  situated  above  the  shelf  break 
instead  of  being  directly  at  the  coast.  The  undercurrent 
has  a  speed  of  ~2  cm/ s  and  has  deepened  to  greater  than 
1500  m  depth. 

Most  obvious  at  day  40  (Eigure  13c)  is  the 
intensification  of  upwelling  south  of  Shark  Bay  and  at 
~29°S.  Although  it  remains  very  localized,  relatively  cold 
water  of  16.5°C  has  been  upwelled  to  the  surface.  The  main 
factor  in  increasing  the  upwelling  is  the  strengthening  of 
the  equatorward  surface  current  to  ~80  cm/s .  Although  the 
current  is  present  along  the  entire  western  and  southern 
coast  of  Australia  it  is  strongest  in  the  northern  end  of 
the  model  domain,  as  expected  due  to  the  stronger  winds  in 
this  region.  As  the  surface  current  increases,  the  core 
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becomes  more  developed  but  continues  to  be  centered  over 
the  shelf  break  (Figure  13d) . 

By  day  60  (Figure  13e) ,  the  southerly  winds  have 
continued  to  enhance  the  equatorward  surface  current,  which 
has  reached  a  maximum  velocity  of  ~120  cm/s .  The  areas  of 
localized  upwelling,  off  Shark  Bay  and  at  29°S,  have  become 
relatively  intense  bringing  ~14°C  water  to  the  surface. 
Temperature  cross-sections  for  the  same  latitudes,  (Figure 
13h  and  13i) ,  respectively,  show  that  the  wide  shelf  to  the 
north  enhances  the  upwelling. 

In  addition  to  being  much  colder  than  the  water 
offshore,  the  surface  elevation  in  the  upwelling  regions  is 
also  lower.  The  alongshore  pressure  gradient  induces  an 
equatorward,  geostrophic  current  which  helps  to  enhance  the 
wind-forced  equatorward  surface  current.  A  relatively 
deep,  poleward  undercurrent  of  ~10  cm/s  is  trapped  by  the 
topography  at  ~1000  m  depth  from  26°S  (Figure  13f)  south  to 
Cape  Leeuwin  (Figure  13g) ,  where  it  becomes  extremely 
narrow . 

Based  on  these  results,  the  addition  of  topography 
traps  and  enhances  the  equatorward  surface  current. 
Simultaneously,  the  topography  hinders  the  development  of 
some  upwelling  regions  (e.g.,  coastal  embayments)  even  with 
upwelling  favorable  winds;  however,  in  areas  where 
upwelling  is  present  (e.g.,  off  capes)  the  upwelling  is 
very  intense.  The  results  of  this  experiment  are 
representative  of  typical  EBC  regions  other  than  the  LCS . 
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F.  THERMOHALINE  GRADIENT  FORCING  OVER  TOPOGRAPHY 

Experiment  6  (see  Table  1)  is  initialized  with  full 
annual  climatology,  but  without  wind  forcing.  It  uses  the 
same  coastline  and  is  run  with  realistic  topography.  It 
introduces  conditions  unique  to  the  LCS . 

By  day  20  (Figure  14a) ,  a  poleward  coastal  current  has 
developed  with  a  maximum  velocity  of  ~90  cm/ s .  The  current 
extends  from  the  northern  edge  of  the  domain  around  the 
southwest  corner  of  Australia  to  the  eastern  edge.  There 
is  also  a  relatively  broad  movement  of  warm,  offshore  water 
to  the  south,  with  temperatures  off  Cape  Leeuwin  reaching 
23°C.  A  typical  cross-section  (Figure  14b)  shows  that  the 
surface  current  is  centered  over  the  shelf  break  and  lies 
inshore  of  an  equatorward  current  which  intensifies  with 
depth . 

By  day  40  (Figure  14c) ,  the  coastal  current  velocity 
has  increased  to  ~140  cm/s .  Warm,  subtropical  water  of  25°C 
has  been  advected  around  Cape  Leeuwin  and  into  the  entrance 
of  the  Great  Australian  Bight  to  the  east  of  the  model 
domain.  Eddies  have  started  to  spin  up  off  the  coast  near 
Cape  Leeuwin  and  are  subsequently  advected  eastward.  An 
eddy  dipole  pair  has  formed  off  Perth.  The  cold  core  half 
of  the  dipole  has  been  created  by  entraining  cold  water.  A 
typical  cross-section  through  a  warm  core  eddy  is  shown  in 
Figure  14d.  Also  showing  is  the  undercurrent,  which 
maintains  its  core  velocity  of  ~20  cm/s  at  ~700  m  depth 
along  the  slope  below  the  poleward  surface  flow. 
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Despite  the  formation  of  many  eddies,  by  day  60 
(Figure  14e) ,  the  Leeuwin  Current  remains  strong  and 
continues  to  advect  warm  water  south  around  Cape  Leeuwin. 
The  current  remains  trapped  by  the  topography  and  is 
centered  above  the  shelf  break  (e.g.,  see  Figure  14f) .  The 
vertical  velocity  shear  caused  by  the  strong,  poleward 
surface  current  and  equatorward  undercurrent  create  large 
meanders  in  the  surface  current,  particularly  in  the  north 
near  Shark  Bay.  As  the  current  travels  around  Cape 
Leeuwin,  it  continues  to  meander  and  create  numerous  eddies 
and  dipole  pairs,  which  are  then  advected  eastward.  One 
dipole  pair,  centered  at  ~35°S,  112°E,  has  separated  from 
the  Leeuwin  Current  and  is  no  longer  under  the  influence  of 
topography.  This  pair  will  propagate  westward  due  to  the  h 
effect.  As  eddies  form,  their  surface  features  deepen  in 
time  and  subsequently  merge  with  subsurface  features  (e.g., 
see  Figure  14f ) . 

Based  on  these  results,  the  addition  of  bottom 
topography  in  this  experiment  isolates  and  traps  the  unique 
thermohaline  gradient  driven,  poleward  current  all  along 
the  coast.  The  topography  aids  the  development  of 
realistic  eddies  and  dipole  pairs,  but  makes  the  core  of 
the  undercurrent  too  deep.  Overall,  the  results  of  this 
experiment  resemble  the  distinct  dynamic  features  that  set 
the  LCS  apart  from  classical  EBCs . 

G.  WIND  AND  THERMOHALINE  GRADIENT  FORCING  OVER  TOPOGRAPHY 

The  final  experiment.  Experiment  7  (see  Table  1), 
includes  all  the  forcing  mechanisms  and  most  closely 
resembles  the  LCS  in  the  annual  sense.  The  model  was 
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initialized  with  full  annual  temperature  and  salinity 
climatology  and  with  a  realistic  bottom  topography  and 
coastline.  In  addition  it  was  forced  with  the  averaged 
annual  wind  field. 

By  day  20  (Figure  15a) ,  a  narrow  poleward  current  has 
formed  with  a  maximum  speed  of  ~90  cm/ s .  This  strong 
surface  current  is  present  all  along  the  shelf  break  and 
advects  warm  subtropical  water  south,  with  temperatures  off 
Cape  Leeuwin  reaching  23°C.  A  cross-section  at  Cape  Leeuwin 
(Figure  15b)  shows  the  core  of  the  current  located  over  the 
shelf  break  with  an  equatorward  undercurrent  of  ~15  cm/s 
centered  near  700  m  depth.  In  regions  where  the 
continental  shelf  is  wider,  such  as  near  Shark  Bay,  the 
main  surface  current  is  formed  as  far  as  50-60  km  offshore. 
This  allows  the  wind  forcing  to  dominate  the  flow  inshore 
of  the  main  current,  generating  an  equatorward  surface 
current  directly  at  the  coast.  A  cross-section  in  such  a 
region  (Figure  15c)  shows  a  poleward  surface  flow  offshore 
above  the  shelf  break,  with  an  equatorward  surface  flow 
inshore  over  the  shallow  continental  shelf. 

By  day  40  (Figure  15d) ,  the  poleward  surface  current 
has  velocity  to  speeds  of  -120  cm/s  and  has  advected 
relatively  warm  24°C  water  south  around  Cape  Leeuwin.  As 
the  surface  velocity  has  increased,  the  meandering  of  the 
current  also  increases  due  to  the  baroclinic  instability 
generated  by  the  increasing  the  vertical  velocity  shear. 
Eddies  are  generated  off  Cape  Leeuwin  and  subsequently 
propagate  eastward  along  the  southern  coast.  The 
equatorward  wind  forcing  pushes  the  cold  subpolar  waters 
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northward,  which  get  entrained  by  the  current  and  help  to 
develop  eddy  dipole  pairs. 

By  day  60  (Figure  15e) ,  the  strength  of  the  current  is 
maintained,  while  it  continues  to  advect  warm,  fresh  water 
south.  Figure  15f  shows  a  surface  plot  of  salinity, 
demonstrating  the  advection  of  fresh  water  by  the  surface 
current  along  the  western  coast.  Many  more  meanders  and 
eddies  have  formed  in  the  current  off  the  west  coast.  The 
addition  of  wind  forcing,  which  opposes  the  poleward 
thermohaline  induced  current,  helps  to  generate  more  well 
defined  eddies.  Two  typical  cross-sections  at  34.3°S 
(Figure  15g)  and  28°S  (Figure  15h)  show  that,  inshore  of  the 
warm  core  eddy,  the  surface  current  hugs  the  shelf  break 
all  along  the  coast.  At  28°S,  note  that  the  wind  forcing 
has  also  created  an  equatorward  surface  current,  with  a 
velocity  of  ~20  cm/s,  inshore  of  the  surface  poleward 
current.  This  equatorward  flow  subsequently  leads  to 
localized  areas  of  upwelling  distinguished  by  slightly 
cooler  waters  along  the  coast,  as  seen  by  the  rising  of  the 
temperature  contours  (e.g..  Figure  151)  and  increasing 
salinity  (e.g., Figure  15f)  over  the  shelf. 

The  results  of  this  final  experiment  show  relatively 
realistic  features  of  the  LCS  including  a  poleward  surface 
current,  equatorward  undercurrent,  meanders  and  eddies,  as 
well  as  regions  of  localized  upwelling.  Unlike  classical 
EBCs  that  are  driven  by  wind  forcing,  the  LCS  is  dominated 
by  the  effects  of  thermohaline  gradients,  which  generate  a 
poleward  (equatorward)  surface  current  (undercurrent)  off 
the  west  coast.  The  thermohaline  gradient  and  bottom 
topography  trap  the  warm  poleward  flow  along  the  coast. 


despite  the  opposing  wind  forcing.  The  wind  forcing 
influences  the  development  of  better  defined  eddies,  and 
generates  equatorward  flow  next  to  the  coast  in  regions 
where  the  poleward  surface  flow  is  farther  offshore. 

H.  DISCUSSION 

Since  Experiment  7  includes  all  the  forcing  mechanisms 
of  the  LCS,  it  is  useful  to  qualitatively  compare  the  final 
results  of  the  model  simulation  with  observational  data  to 
determine  if  the  model  developed  a  current  system 
representative  of  the  LCS.  Because  the  study  is  not  based 
on  actual  data  but  instead  is  a  process-oriented  study, 
direct  comparisons  to  observational  data  cannot  be  made; 
however  the  dynamics  of  the  model  can  be  investigated  to 
determine  if  there  are  qualitatively  similar  to 
observations  of  the  LCS. 

The  results  from  Experiment  7,  day  60  (Eigure  15e)  , 
will  be  used  in  the  comparison.  The  general  features  of 
the  Leeuwin  Current  apparent  in  the  satellite  images 
(figures  16  through  18)  are  similar  to  the  results  found  in 
the  final  model  experiment.  The  strong  thermohaline 
gradient  creates  an  onshore  geostrophic  flow,  which  then 
forms  a  narrow  poleward  current  advecting  warm  subtropical 
waters  southward  around  Cape  Leeuwin  and  into  the  Great 
Australian  Bight.  All  three  satellite  images  show  evidence 
of  meanders,  jets  and  eddies  with  varying  length  scales 
which  are  consistent  with  the  model  results.  Studies  have 
shown  that  the  LCS  has  large  spatial  variations  due  to 
intensification  of  instabilities  (Batteen  and  Butler,  1998; 
Cresswell,  1981;  and  Pearce  and  Giffiths,  1991).  Current 
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velocity  speeds  of  ~120  cm/s  at  day  60  of  Experiment  7, 
agree  with  surface  current  measurements  greater  than  100 
cm/s  made  by  the  R/V  Franklin  at  Cape  Leeuwin  (Cresswell 
and  Peterson,  1993)  .  Overall  the  model  results  for  this 
study  are  supported  by  their  resemblence  to  both  satellite 
images  and  previous  studies  of  the  LCS . 

A  final  comparison  of  all  the  experiments  is  made  to 
determine  the  role  of  each  forcing  mechanism  in  generating 
the  LCS.  For  this  comparison  the  last  day  of  each 
experiment (excluding  Experiment  4,  the  pressure  gradient 
force  error),  i.e.  day  60  (Figures  9e,  lOe,  lie,  13e,  14e 
and  15e) ,  will  be  analyzed.  The  obvious  dominant  forcing 
mechanism  for  the  LCS  is  the  thermohaline  gradient  which 
has  a  strong  poleward  component.  The  gradient  sets  up  an 
onshore  geostrophic  flow,  which  when  it  approaches  the 
coast,  is  forced  poleward  into  a  strong,  narrow  surface 
poleward  current.  The  flow  is  steered  down  the  coast  and 
as  it  reaches  Cape  Leeuwin  it  forms  meanders  and  eddies. 
In  regions  where  the  Leeuwin  Current  leaves  the  coast,  the 
wind  forcing  can  control  the  flow  and  create  an  equatorward 
flow,  which  is  upwelling  favorable.  Since  the  wind  opposes 
the  main  surface  poleward  current,  the  southerly  wind 
generates  considerable  shear,  which  aids  in  the  development 
of  better  defined  eddies. 

While  wind  and  thermohaline  forcing  are  the  main 
driving  mechanisms  for  the  current,  without  topography,  the 
results  from  the  flat  bottom  experiments  show  broad 
currents,  which  do  not  accurately  represent  the  LCS.  In 
terms  of  current  width,  the  topographic  beta  effect  has 
greater  influence,  by  trapping  the  current  along  the  coast. 
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than  the  planetary  beta  effect,  which  broadens  the  current. 
In  the  flat  bottom  case,  the  current  is  steered  by  the 
coastline.  When  topography  is  added,  the  current  is 
steered  by  the  shelf  break,  frequently  leaving  the  coast  as 
it  follows  the  200m  depth  contour  southward  to  Cape  Leeuwin 
and  eastward  into  the  Great  Australian  Bight.  A  comparison 
of  the  wind  only  experiments  shows  that  topography  reduces 
the  upwelling  that  was  present  along  the  entire  coast  in 
the  flat-bottom  experiment  and  enhances  localized  areas 
just  north  of  Cape  Leeuwin  and  of  30°S.  When  the  current 
is  not  trapped  along  the  coast,  such  as  in  the  north  near 
Shark  Bay  and  just  north  of  Cape  Leeuwin,  the  wind  forcing 
generates  an  equatorward  current  which  consequently 
produces  upwelling.  The  presence  of  upwelling  in  these 
areas  is  important  to  the  biological  productivity  of  the 
region,  which  is  important  to  the  fishing  and  tourism 
industries . 
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IV.  SUMMARY 


The  objective  of  this  process -oriented  study  was  to 
build  on  previous  studies  and  to  investigate  the  roles  of 
wind  forcing,  thermohaline  gradients,  and  bottom  topography 
in  the  Leeuwin  Current  System.  The  results  of  several 
numerical  experiments  (see  Table  1)  using  the  POM,  a 
bottom-following  sigma  coordinate  model,  was  explored.  The 
POM  was  chosen  for  this  study  because  it  has  been  widely 
used  to  simulate  coastal  processes  associated  with 
continental  shelf  flows  and  bottom  boundary  layer  dynamics. 
In  all  experiments  a  beta-plane  was  used.  Experiment  1 
(Experiment  2)  investigated  the  effect  of  annual  wind 
forcing  (full,  annual  climatology)  over  a  flat  bottom. 
Experiment  3  was  run  with  both  annual  wind  forcing  and 
full,  annual  climatology  over  a  flat  bottom.  The  results 
of  Experiment  4  were  used  to  determine  the  pressure 
gradient  force  error  created  by  the  addition  of  bottom 
topography  to  the  model.  Experiment  5  (Experiment  6) 
investigated  the  role  of  annual  wind  forcing  (full,  annual 
climatology)  over  topography.  The  final  experiment. 
Experiment  7,  was  run  with  full,  annual  climatology  and 
wind  forcing  over  topography. 

In  Experiment  1  the  annual  wind  forcing,  without  a 
thermohaline  gradient,  forced  a  weak  equatorward  surface 
current.  It  also  produced  classical  EEC  features  including 
upwelling,  meanders,  eddies  and  a  poleward  undercurrent. 
The  upwelling  was  evident  along  the  entire  western  coast  of 
Australia,  but  was  most  intense  near  Shark  Bay. 
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Experiment  2  removed  the  wind  forcing  but  introduced 
the  annual  temperature  and  salinity  gradients,  over  a  flat 
bottom.  The  effects  of  the  thermohaline  gradient  developed 
a  strong,  broad  poleward  surface  current  more 
characteristic  of  the  LCS .  Planetary  beta  effects  created 
westward  propagation  of  the  current  increasing  its  offshore 
extent.  Warm,  subtropical  waters  were  advected  southward 
around  Cape  Leeuwin.  This  experiment  also  generated 
meanders  and  eddies  south  of  Perth,  as  well  as  an 

equatorward  undercurrent. 

The  final  flat  bottom  experiment.  Experiment  3, 

combined  the  effects  of  annual  wind  forcing  and  the 
thermohaline  gradient.  The  results  resembled  those  of 
Experiment  2,  showing  that  the  thermohaline  gradient  is  the 
dominating  forcing  mechanism  in  the  LCS.  The  role  of  the 
wind  was  to  slow  the  speed  of  the  poleward  current  and 

enhance  the  generation  of  eddies.  The  wind  also  created 
localized  upwelling  regions  along  the  coast  where  the 
poleward  flow  left  the  coastline.  Einally,  the  wind 
forcing  advected  subpolar  water  northward,  which  frequently 
became  entrained  by  current  meanders  and  subsequently 
formed  the  cold  core  half  of  eddy  dipoles. 

Experiment  4  introduced  topography  but  was  run  with  no 
wind  or  thermohaline  gradient  forcing  to  determine  the 

pressure  gradient  force  error.  It  showed  that  the  velocity 
errors  inherent  in  a  sigma  coordinate  model  could  be 
reduced  from  ~1  m/s  to  ~.01  m/s  using  three  techniques. 
The  techniques  used  were:  1)  smoothing  the  topography,  2) 
using  the  highest  possible  resolution,  particularly  along 
the  shelf  break,  and  3)  subtracting  the  area-averaged 
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density  before  computation  of  the  baroclinic  integral.  The 
results  showed  that  the  highest  velocities  induced  by  the 
pressure  gradient  force  error  were  directly  related  to 
regions  were  the  slope  parameter  was  the  largest,  i.e.  the 
shelf  slope  and  coastline. 

Experiment  5  explored  the  role  of  annual  wind  forcing, 
without  annual  climatology,  over  bottom  topography.  The 
topography  had  the  effect  of  trapping  the  equatorward 
current  along  the  coast,  increasing  its  maximum  velocity 
when  compared  to  the  flat  bottom  case  (Experiment  1) . 
Although  the  topography  limited  the  extent  of  upwelling  to 
regions  north  of  30°S,  the  intensity  of  the  upwelling  was 
increased.  Relatively  cool  water  of  14.5°C  was  upwelled 
near  Shark  Bay. 

In  Experiment  6,  the  thermohaline  gradient  forcing  was 
isolated  over  bottom  topography.  Similar  to  Experiment  5, 
the  addition  of  topography  opposed  the  planetary  beta 
effect  and  trapped  the  current  along  the  coast  creating  a 
strong,  narrow  poleward  current.  The  shelf  break  steered 
the  current  southward  along  the  coast  and  around  Cape 
Leeuwin  into  the  Great  Australian  Bight  region.  The 
trapped  coastal  current  had  greater  instabilities  than  in 
the  flat  bottom  experiments  resulting  in  more  meanders, 
offshoots  and  eddies. 

Experiment  7  provided  the  best  representation  of  the 
LCS  because  it  included  annual  wind,  annual  climatology, 
and  bottom  topography.  In  terms  of  location  and  velocity, 
it  generated  the  most  realistic  equatorward  surface  current 
and  poleward  undercurrent.  The  results  of  this  experiment 
show  that,  in  an  annual  sense,  the  effects  of  thermohaline 
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gradients  dominate  and  set  the  system  apart  from  classical 
EBCs.  The  major  role  of  the  wind  is  to  slow  the  poleward 
surface  flow,  enhance  eddy  spin-up  and  create  localized 
upwelling  regions.  Bottom  topography  is  shown  to  play  an 
important  role  in  intensifying  and  trapping  currents  near 
the  coast,  in  weakening  the  subsurface  current  and  in 
intensifying  eddies  off  capes.  Overall,  Experiment  7,  the 
most  complete  experiment  of  the  study,  compared  well  to 
observations  and  provided  the  most  accurate  representation 
of  the  LCS .  Although,  in  the  annual  sense,  the  wind  is  an 
insignificant  factor  compared  to  the  thermohaline  gradient, 
it  is  the  seasonality  of  the  wind  that  determines  when  the 
Leeuwin  Current  flows.  further  research  using  seasonal 
winds  and  seasonal  thermohaline  forcing  is  recommended  due 
to  the  seasonal  nature  of  the  observed  Leeuwin  Current. 
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Figure  1 .  The  model  domain  for  the  Leeuwin  Current 
System  (LCS)  is  bounded  by  39°S  to  24.5°S,  109°E  to 

121°E.  The  model  domain  has  a  closed  boundary  along 
the  entire  coast  and  four  open  boundaries. 
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Figure  2a.  Original  topography  (meters)  in  3-D  (from 
Sandwell  and  Smith,  1996)  with  a  resolution  of  2 
minutes  (1/30  degree).  Contours  depict  the  shoreline, 
200m  and  2500m  isobaths. 
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Figure  2b.  Original 
Sandwell  and  Smith, 
minutes  (1/30  degree) 


topography  (meters)  in  2-D 
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Figure  3.  Resolution  grid  lines  with  every  fifth  grid 
line  plotted.  (10  by  10  km  res.  with  two  bands  of  10 
by  3  km  and  3  by  10  km  resolution) . 
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Smoolhed  Topograpy 


Figure  4a.  Smoothed  topography  (meters)  obtained 
after  applying  a  two-dimensional  Gaussian  filter  and 
reassigning  depths  greater  than  2500  m  to  2500  m. 
Contours  depict  the  shoreline,  200  m  and  2500  m 
isobaths . 
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Smoothed  Topography 
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Figure  4b.  Smoothed  topography  (meters)  in  2-D 

obtained  after  applying  a  two-dimensional  Gaussian 
filter  and  reassigning  depths  greater  than  2500  m  to 
2500  m. 
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Figure  5.  Plot  of  the  21  sigma  levels. 
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Figure  6.  Annual  climatological  surface  temperature 
(°C)  obtained  from  Levitus  and  Boyer  (1994) . 
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Figure  7.  Annual  climatological  surface  salinity  from 
Levitus  et  al .  (1994). 
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Figure  8.  Annual  average  wind  in  m/s 

climatological  ECMWF  winds  obtained  from  Trenbert 
al..  1990. 


Figure  9a.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  1  on  day  20. 
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Figure  9b.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  1  day  20. 
Blue  is  equatorward  (north) . 
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Figure  9c.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  1  on  day  40. 
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Figure  9d.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  1  day  40. 
Blue  is  equatorward  (north) . 
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Figure  9e .  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  1  on  day  60. 
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Figure  9f.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  1  day  60. 
Blue  is  equatorward  (north) . 
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Figure  10a.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  2  on  day  20. 
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Figure  10b.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  2  day  20. 
Blue  is  equatorward  (north) . 
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Figure  10c.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  2  on  day  40. 
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Figure  lOd.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  2  day  40. 
Blue  is  equatorward  (north) . 
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Figure  lOe.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  2  on  day  60. 


53 


Velocity 


Longitude 


Figure  lOf.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  2  day  60. 
Blue  is  equatorward  (north) . 
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Figure  11a.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  3  on  day  20. 
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Figure  11b.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  3  day  20. 
Blue  is  equatorward  (north) . 
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Figure  11c.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  3  on  day  40. 
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Figure  lid.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  3  day  40. 
Blue  is  equatorward  (north) . 
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Figure  lie.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  3  on  day  60. 
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Figure  Ilf.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  3  day  60. 
Blue  is  equatorward  (north) . 
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Figure  llg.  Cross-section  of  north-south  velocities 
(m/s)  at  Perth  (32.5°S)  for  Experiment  3  day  60.  Blue 
is  equatorward  (north) . 
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Figure  12.  Surface  velocity  error  (m/s)  due  to  the 
pressure  gradient  force  error  on  day  60  for  Experiment 

4  . 
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Figure  13a.  Surface  temperature 
vectors  for  Experiment  5  on  day 
subsequent  horizontal  plots,  the 
shown  by  the  solid  line. 


(°C)  and  velocity 
2  0.  In  this  and 
200  m  isobath  is 


63 


Velocity 


Longitude 


Figure  13b.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  5  day  20. 
Blue  is  equatorward  (north) . 
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Figure  13c.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  5  on  day  40. 
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Figure  13d.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  5  day  40. 
Blue  is  equatorward  (north) . 
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Figure  13e.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  5  on  day  60 . 
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Figure  13f.  Cross-section  of 
(m/s)  at  26°S  for  Experiment 
equatorward  (north)  . 


north-south  velocities 
5  day  60.  Blue  is 
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Figure  13g.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  5  day  60. 
Blue  is  equatorward  (north) . 
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Figure  13h.  Cross-section  of  temperature  (°C)  at  32°S 
for  Experiment  5  day  60,  with  a  contour  intervals  of 
two  degrees. 
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Figure  13i.  Cross-section  of  temperature  (°C)  at  26°S 
for  Experiment  5  day  60,  with  a  contour  intervals  of 
two  degrees. 
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Figure  14a.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  6  on  day  20. 
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Figure  14b.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  6  day  20. 
Blue  is  equatorward  (north) . 
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Figure  14c.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  6  on  day  40. 
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Figure  14d.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  6  day  40. 
Blue  is  equatorward  (north) . 
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Figure  14e.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  6  on  day  60. 
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Figure  14f.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  6  day  60. 
Blue  is  equatorward  (north) . 


77 


Figure  15a.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  7  on  day  20. 
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Figure  15b.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  7  day  20. 
Blue  is  equatorward  (north) . 
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Figure  15c.  Cross-section  of  north-south  velocities 
(m/s)  at  27°S  for  Experiment  7  day  20.  Blue  is 
equatorward  (north) . 
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Figure  15d.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  7  on  day  40. 
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Figure  15e.  Surface  temperature  (°C)  and  velocity 
vectors  for  Experiment  7  on  day  60. 
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Figure  15f.  Surface  salinity  and  velocity  vectors  for 
Experiment  7  on  day  60 . 
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Figure  15g.  Cross-section  of  north-south  velocities 
(m/s)  at  Cape  Leeuwin  (34.3°S)  for  Experiment  7  day  60. 
Blue  is  equatorward  (north) . 


84 


Velocity 


Longitude 


Figure  15h.  Cross-section  of  north-south  velocities 
(m/s)  at  28°S  for  Experiment  7  day  60.  Blue  is 
equatorward  (north) . 
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Figure  15i.  Cross-section  of  temperature 
for  Experiment  7  day  60,  with  a  contour 
two  degrees. 
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Figure  16. 
temperatures 
showing  the 
(red/ orange) 
green/blue . 


Satellite  image  of  surface  water 
off  Western  Australia  in  August  1987, 
warm  waters  of  the  Leeuwin  Current 
and  the  cooler  offshore  water  in 
The  white  and  mottled  blue  areas  are 


clouds,  and  the  black  line  represents  the  edge  of  the 
continental  shelf (CSIRO,  Marine  Research) . 
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Figure  17.  Satellite  image  of  surface  water 
temperatures  off  Western  Australia  in  September  1994, 
showing  the  warm  waters  of  the  Leeuwin  Current 
(red/orange)  and  the  cooler  offshore  water  in 
green/blue.  The  white  and  mottled  blue  areas  are 
clouds,  and  the  black  line  represents  the  edge  of  the 
continental  shelf (CSIRO,  Marine  Research) . 
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Figure  1 8  . 
temperatures 
showing  the 
(red/ orange) 
green/blue . 

and  the  black  line  represents  the  edge 
continental  shelf (CSIRO,  Marine  Research) . 


Satellite  image  of  surface  water 
off  Western  Australia  in  August  1987, 
warm  waters  of  the  Leeuwin  Current 
and  the  cooler  offshore  water  in 
The  white  and  mottled  areas  are  clouds, 

of  the 
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Exp  # 

Annual  Wind 

Climatology 

Topography 

1 

YES 

HORIZONTALLY 

AVERAGED 

NO 

2 

NO 

EULL 

NO 

3 

YES 

EULL 

NO 

4 

NO 

HORIZONTALLY 

AVERAGED 

YES 

5 

YES 

HORIZONTALLY 

AVERAGED 

YES 

6 

NO 

EULL 

YES 

7 

YES 

EULL 

YES 

Table  1.  Summary  of  experimental  design. 
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Level 

Depth  (m) 

Level 

Depth  (m) 

Level 

Depth  (m) 

1 

0 

12 

300 

23 

1400 

2 

10 

13 

400 

24 

1500 

3 

20 

14 

500 

25 

1750 

4 

30 

15 

600 

26 

2000 

5 

50 

16 

700 

27 

2500 

6 

75 

17 

800 

28 

3000 

7 

100 

18 

900 

29 

3500 

8 

125 

19 

1000 

30 

4000 

9 

150 

20 

1100 

31 

4500 

10 

200 

21 

1200 

32 

5000 

11 

250 

22 

1300 

33 

5500 

Table  2.  Vertical  levels  and  depths  used  by  Levitus 
and  Boyer  (1994)  and  Levitus  et  al .  (1994) 
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Table  3.  Values  of  sigma  levels 
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